S ; Criticality and Condensation in a Non-conserving 
o . Zero Ranee Process 

CN ■ ° 

^ ; A. G. M. R. Evans^, E. Levine^ and D. Mukamel'' 

I ^Center for Stochastic Processes in Science and Engineering, Department of Physics, 

00 . Virginia Tech, Blacksburg, VA 24061-0435, USA. 

^SUPA, School of Physics, University of Edinburgh, 

; Mayfield Road, Edinburgh EH9 3JZ, UK. 

•^Center for Theoretical Biological Physics, 

g I University of California at San-Diego, La JoUa, CA 92093, USA. 

''Department of Physics of Complex Systems, Weizmann Institute of Science, 

, Rehovot, Israel 76100. 
> . 

' E-mail: aangelSowl.phys.vt.edu.m. evaiis9ed.ac.uk, levine@ctbp.ucsd.edu, 

, david.mukcmieKSweiznicLnn.ac. il 

S: 

^ ' Abstract. The Zero-Range Process, in which particles hop between sites on a 

jH . lattice under conserving dynamics, is a prototypical model for studying real-space 

O I condensation. Within this model the system is critical only at the transition point. 

Here we consider a non-conserving Zero-Range Process which is shown to exhibit 
generic critical phases which exist in a range of creation and annihilation parameters. 
The model also exhibits phases characterised by mcsocondensates each of which 
contains a subextensive number of particles. A detailed phase diagram, delineating 
\ the various phases, is derived. 
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1. Introduction 

There exist many systems exhibiting real-space condensation under nonequihbrium 
steady state conditions. Examples include jamming in traffic flow, granular clustering, 
wealth condensation and gelation in networks [1]. Real-space condensation implies that 
a finite fraction of some conserved quantity, for example density, condenses onto a single 
lattice site or a small region in space. In general one is interested in the distribution of 
the conserved quantity and the emergence of the condensate as some parameter, often 
the density, is varied. 

The Zero-Range Process (ZRP) [2] is a generic model which exhibits condensation 
in its nonequihbrium steady state [3l H] . Moreover it has the convenient property that its 
steady state is known and has a simple factorised form. Thus the ZRP provides a simple 
exactly solvable model within which general features of condensation may be studied 
and analysed [1]. For example, recent developments include: extensions of the model to 
several conserved quantities [HlEllTj, open boundary conditions [8], and sitewise disorder 
[9l[T0l[TT]; the study of current fiuctuations [12] and traffic modelling [T3| [T^. 

The ZRP is usually defined on a lattice of L sites where each site may be occupied 
by any integer number of particles. The dynamics is defined by hopping rates u{n) with 
which a particle hops from a site occupied by n particles. In various versions of the 
dynamics the particle may hop to different allowed destination sites. For example on 
a regular lattice particles hop to nearest neighbour sites whereas on a fully-connected 
geometry a particle can hop to any other site with equal probability. Clearly, the total 
number of particles, A^, is conserved under the ZRP dynamics. Condensation occurs 
when in the large A^, L limit, with the density p = N/L held fixed, a finite fraction of 
the particles condenses onto a single lattice site. More recently, an example where the 
condensate has a non-zero spatial extent has been studied |15j . 

Condensation is revealed in the single-site occupation probability distribution, 
p{n). A characteristic case is when the hopping rate has the asymptotic, large n form 
u{n) ^ 1 + b/n. For b > 2, the model exhibits a condensation transition at a critical 
density Pc- For p < pc, p{n) has the form 



where /i is positive and is a function of the density. Since there is a characteristic 
occupation this is referred to as the fiuid phase. 

As the density increases towards the critical value pc, fi tends to zero and the 
resulting distribution becomes a power law at p = pc- For p > pc an extra piece of 
p(n), representing a single condensate, emerges centered at n = L(p — pc) [16]. Thus in 
the condensed phase a critical fiuid co-exists with a condensate containing the excess 
density. A pure power-law distribution only holds at the critical density which is a 
common feature in many systems exhibiting phase transitions. For b < 2 there is no 
condensation transition since any density can be achieved by choosing p to be suitably 
small in ([1]). 




(1) 
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In this work we consider more general dynamical processes which could allow for 
more complex condensation phenomena. Examples of the phenomena we have in mind 
are the existence of criticality in an extended region of the phase diagram rather than 
at an isolated point (typically referred to as self-organised criticality), and condensation 
into a large number of condensates each containing a subextensive number of particles. 
In our generalised dynamics we introduce non- conserving processes with creation and 
annihilation rates. One mechanism to suppress a single extensive condensate is for the 
annihilation of particles to occur preferentially at highly occupied sites. We find that 
an appropriate choice of creation and annihilation rates leads to a whole critical region 
of the phase diagram where the occupation distribution decays algebraically at large 
occupations. Since the density is not conserved the phase diagram is given in terms of 
the parameters of the creation and annihilation rates. 

Let us briefly summarise the phases which are exhibited (see Fig. [1]). As might be 
expected, for imbalanced creation and annihilation rates we find regimes with vanishing 
or diverging density. These two regimes are separated by phases where the observed 
density is equal to what would be the critical density on a conserving system and 
power-law occupation distributions are exhibited. One of these phases. Critical Phase 
A, exhibits a pure power-law distribution for p{n) with a cut-off diverging with system 
size. However another region in parameter space exists where the distribution exhibits, 
in addition to the power-law decay, a broad and weak peak at high occupations. The 
height of the peak scales algebraically with system size, rather than exponentially, as 
does the width. This peak may correspond to many 'mesocondensates' each of which 
contains a subextensive number of particles. The number of these mesocondensates is 
also subextensive thus they occupy a vanishing fraction of the sites. We shall give a 
detailed analysis of this weak peak and how it leads to two phases: Critical Phase B 
where the weak peak does not contribute to the global density which is pc and the Weak 
High Density Phase where the contribution of the weak peak to the global density is 
dominant. 

In a recent publication [l7j we gave a brief account of the non- conserving ZRP and 
its relation to self-organised criticality and the dynamics of rewiring networks. Here we 
focus on the properties of the non-conserving ZRP and present a detailed analysis of 
the phase diagram. 

The structure of the paper is as follows: in Section [2] we define the non-conserving 
ZRP that we study; in Section [3] we present a detailed analysis of the phase diagram 
via a mean-field approximation and we compare the results to numerical simulations on 
fully-connected and one-dimensional lattices; we conclude with a discussion in Section HI 

2. Definition of non-conserving ZRP 

We now define the non-conserving zero-range process which we study in this work. The 
lattice contains L sites labelled / = 1, . . . , L. The number of particles at site / is ri;. The 
dynamics are defined by the following three processes: particles hop from a site with n 
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particles with rate 

u{n) = (l + ^) , (2) 

where the step function 6{n) is defined as 

/ 1 for n > , . 

'^") = io for n = 0, 

particles are created at a site with rate 

and particles evaporate from a site containing n particles with rate 

a{n) = [^^y e{n) . (5) 

In these rates the indices k and s are positive. The creation rate at a site (jlj) provides a 
weak drive which decreases with system size. The particle annihilation rate at a site ([5]) 
increases with the number of particles at the site and provides a mechanism by which 
an extensive condensate may be suppressed. 

As noted in the introduction, condensation in the conserving ZRP occurs when 
b > 2 and in the following we focus primarily on this range of b. For the specific choice 
of hopping ([2]) the critical density is known to be [181 [H] 



Note that in the non-conserving case and consequently the density p fiuctuate in 
time. 

The precise nature of the lattice and the definition of the sites to which particles are 
allowed to hop from a given site do not affect the qualitative picture emerging from this 
study (as long as the lattice is homogeneous with all sites having the same hop rates). 
To be specific we consider a fully-connected lattice where a particle can hop to any other 
lattice site. This is a convenient choice, since on the fully-connected lattice the mean- 
field approach with which we treat the model analytically is expected to become exact 
in the large system limit. In addition, we shall also present numerical data from one 
dimensional systems with totally asymmetric nearest neighbour hopping which indicate 
that the results of the mean-field analysis are relevant to this case as well. 

In a Monte-Carlo simulation the dynamics defined in ([2H5]) are conveniently 
implemented as follows: At each update 

(i) Select a site at random. 

(ii) Generate a random number uniformly distributed between zero and the sum of the 
maximum possible rates of the three processes, i.e., c + max(a(n)) -|- max(M(?T,)). 

(iii) If the random number falls in the range (0, c) create a particle at the site. 

(iv) If the random number falls in the range (c, c -|- a(n)), where n is the number of 
particles on the site, evaporate a particle from the site. 
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(v) If the random number falls in the range (c + max(a(n)),c + max(a(n)) + u{n)) 
remove a particle from this site and place it at another randomly selected site. 

Note, as the evaporation rate a{n) can diverge, a cut-off max(a(r2)) must be imposed 
artificially and chosen large enough that, in practice, the dynamics is not affected. 

The steady state of the model is fully described by the probability distribution 
P{ni, 712, ■ ■ ■ , ^l) over all possible configurations. In contrast to the conserving ZRP [1], 
the steady-state distribution of the model ([2H5]) does not factorize generally. However, 
for the case of a fully-connected lattice we expect factorisation to take effect in the 
thermodynamic limit L oo. In the following analysis we apply a mean field 
approximation in which the steady state distribution is replaced by a factorized form, 
i.e. P(ni, 77,2, . . . , n^) ]\i=iP{ni). For the fully-connected lattice we expect this 
approximation to become exact in the thermodynamic limit. 

Within this approximation, the master equation is given by 

= [u{n + 1) + a{n + 1)] p{n + 1) - [A + c] p{n) (7) 
— {[u{n) + a{n)] p{n) — [A + c] p{n — 1)} 6{n) . 

where p{n) is the probability that a site contains n particles. The step function 9{n) 
ensures that ([7j) holds for all n > 0. The 'current' A into a site due to the hopping 
process is given by 

oo 

A = ^ u{n)p{n) . (8) 

n=l 

To understand ([7]) note that u{n) + a{n) is the total rate at which a site with n particles 
loses a particle and A + c is the total rate at which a site gains a particle. 
In the steady state ([7]) becomes 

= [u{n + 1) + a{n + 1)] p{n + 1) - [A + c] p{n) (9) 

— {[u{n) + a{n)] p{n) — [A + c] p{n — 1)} 6{n). 

which may be iterated to obtain the steady-state p{n) 

p{n) = ,^-J:^±^L^p(0) . (10) 

To fix p(0) and A, p{n) must satisfy the following constraints of normalisation and 
creation-annihilation balance respectively 

oo 

EpW =1 (11) 

n=0 
oo 

E a{n)p{n) = c . (12) 

n=l 

Only two of equations ([8]), (ITT]) and (fT2|) are independent, as can be seen by noting 

[a{n) + u{n)] p{n) = (A + c)p{n — 1) , (13) 
which, when summed, becomes 

oo oo 

[a{n) + u{n)] p{n) = (A + c) ^ p{n) = A + c . (14) 

n=l n=0 
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3. Analysis of phase diagram of non-conserving ZRP 

We now determine the asymptotic, large L behaviour of A (and therefore p{n)) required 
to satisfy ([8]) and f|T2|) . We identify a total of five phases in the s-k plane: low density 
phase, strong and weak high density phases and a critical region which itself consists of 
two distinct critical phases as we shall explain below. 

4 
3 

1 



12 3 4 

k 

Figure 1. Typical phase diagram for the non-conserving ZRP model, shown in the 
k-s plane for 6 = 3. The parameters k,s,b are defined in ([2H5]). The labelling of the 
phases corresponds to Low Density Phase (LD), Critical Phase A (CA), Critical Phase 
B (CB), Weak High Density Phase (WHD) and Strong High Density Phase (SHD). 




Identifying the five phases is most conveniently done by inspecting the balance 
equation (fT2l) which becomes on inserting the expressions (fTOl) for p{n), (jlj) for c and 
([5]) for a 

oo 
n=l 

oo 

= p(0) exp 

n=l 

3.1. Low density phase: s > k 

For s > k the LHS of (ITSl) tends to zero as the system size L tends to infinity. This 
requires that p{n) is a rapidly decreasing function of n. To satisfy the balance equation 
f fTSjl A + c has to be small. Then to leading order in L, A + c ^ A ~ L'^"'* 

p{n) ~ L-^^-'^)" . (16) 

Thus the density p ~ p{l) ~ V''^ which tends to zero as the system size tends to oo. 
The mean total number of particles increases sublinearly with system size as 
and so we are in a low density phase when s > k. 



n ln(A + c) 



n 

E 

m=l 



ln(a(m) + u{m)) 



(15) 
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3.2. Analysis of s < k region 

In this case, since s < k, the LHS of ( fT5|) diverges with L, so the sum on the RHS must 
also diverge and must be dominated by terms at large n. This can only happen if A + c 
approaches one for large L. Otherwise, if the limiting value of A + c were less than one 
the sum would not diverge, whereas if the limiting value were greater than one the sum 
would diverge even for finite L. Therefore we write 

A + c = 1 + h{L) - e'^(^) , (17) 

where h{L) ^ as L ^ oo. Since we are interested in the large n behaviour of p{n), 
keeping leading order terms in n, we have ln(a(m) + u{m)) c:^ b/m + {m/ L)^ and then 
approximating the sum by an integral gives, for large n. 

In [aim) + M(m)] ^ h\n{n) + 



m=l 



n 

{k + l)L^ ■ 



Therefore, from fllUp . the asymptotic form of p{n) is 

h{L)n 



pH ~ 4: exp 



{k + l)L>' 

The rest of the analysis amounts to determining more precisely the large L behaviour 
of h{L) required to satisfy the balance equation f[T^ in the various regions of the s-k 
plane. 

3.3. Strong High Density Phase: s < k/{k + 1) 

The form of ( |T8l) suggests that if h{L) is positive, p{n) may be sharply peaked at some 
value n* which dominates p{n). By sharply peaked we mean that for e arbitrarily small, 
/"*^/dn p{n) — >• 1 as L — > oo. Then the balance equation (IT^ reduces to L'^~^ ~ (n*)^ 
which implies 

n* ~ L^-'/'' . (19) 

This condition may be used to determine h{L). Maximising the argument of the 
exponential in f|T8l) yields n* = Lh^^^{L). Therefore comparing with f|T9|l we require 

h{L) ~ . (20) 

However, for the distribution to be sharply peaked at n* we require that the argument of 
the exponential in ffTSjl diverges for large L at n* . Thus, for example, n*h{L) ~ j^^-s-s/k 
should diverge which implies 

In this phase, note that A + c = 1 + 0{L *) approaches one from above. Also note that 
the number of particles is super extensive and the mean density is n* which diverges as 
f|T9|) . The sharply peaked distribution p{n) is rather different from that of a conserving 
ZRP at high density, where a critical fluid and a condensate piece coexist. In contrast, 
in the present case all sites contain 0{n*) particles and we refer to this as the Strong 
High Density Phase. 
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3.4- Intermediate Regime 

We now consider the intermediate regime k/{k + 1) < s < k and show that p{n) 
approaches a power-law distribution. A careful analysis shows that, in fact, this region 
may be divided into three phases: two of these are critical in the sense that the density is 
given by pc, the critical density of the conserving ZRP. The first critical phase (Critical 
Phase A) has a power-law distribution p{n) with a cut-off at large n. The second critical 
phase (Critical Phase B) has a power-law distribution together with a broad and weak 
peak at large occupations. In the thermodynamic limit this peak does not contribute to 
the global density and thus the global density is critical. Finally we have a phase similar 
to Critical Phase B but where the broad peak dominates the global density which now 
diverges with L. We refer to this as the Weak High Density Phase. 



3.4-1 ■ Critical Phase A {kb/{k + 1) < s < k) In the case where X + c approaches 
one from below, namely when h{L) in ( |T7I) is negative, the distribution ( ITSj) becomes 
a power law with a large n cut-off at n ~ mm[l / \h(L)\, L'^^^''^^^]. This cut-off goes to 
infinity with L resulting in a power-law distribution in the thermodynamic limit. 

We analyse the region in the s-k plane where such a behaviour is manifested. 
Setting h{L) ~ —gL~^, where g is a. positive constant, and replacing the sum by an 
integral, the balance equation (fT5l) takes the form 



dn ^ 



exp 



n 



n 



(22) 



(A; + 1)^^= 

In the case k > h — 1 the integral is dominated by large n. The cut-off in the large n 
contribution will be given by L^, as long a.s x < k/ {k + 1). The asymptotic behaviour 
of (1221) becomes 



Thus, X is given by 



k 



X 



k-b + 1 

To be consistent with the requirement x < k/{k 
satisfy 

hk 



1) one finds from ([2 



s > 



k 



(23) 
(24) 

that s has to 
(25) 



When the equality holds, both terms in the exponential in (1221) are relevant. To 
summarise, in the range kb/ (k + l) < s < k the asymptotic behaviour of p{n) is a power 
law with an exponential cut-off, the cut-off point tending to infinity as the system size 
tends to infinity: 

n 



p{n) 



~ n exp 



-9 



(26) 



where x is given by (12^ . and g is an undetermined positive constant. For this phase to 
exist we require k > b — 1. 
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3.4-2. The region k/{k + 1) < s < kh/{k + 1) This still leaves the behaviour in the 
region A;/(A; + 1) < s < kh/ {k + 1) to be determined. We expect that some kind of local 
maximum in p[n) is needed in order to cope with the increased creation rate in this 
region compared with Critical Phase A. However, a sharp peak in p{n) is ruled out since 
that would correspond to the Strong High Density phase. Instead we shall show that a 
'weak peak' emerges in p{n) whose height scales algebraically with L as opposed to the 
sharp exponential peak in the high density phase (see Figure [2]). We now demonstrate 
the existence of such a peak in the distribution in the region k/{k + l) < s < kh/ {k + 1), 
and show how it satisfies the balance equation (fT5l) . 



p{n*) - 




'^niin 



Figure 2. Schematic illustration of a weak peak in p{n). Note that both the 
height p{n*) and width An* scale algebraically with system size L but the weight. 
w = An*p{n*) of the peak vanishes. The dip in the distribution to the left of the weak 
peak is denoted n,^i^. 



Let us write the distribution (1181) as 

p{n) ~ e'^^") (27) 



where 



^Pin) = h{L)n - -b\nn. (28) 

We now seek a function h{L) which produces a maximum at n* and for which 
ip{n*) ~ InL. This would result in p{n*) scaling as a power of L. In order for the 
second term in fl28l) to scale as In L we require n* ~ (LMnL)i/(fc+i). Then for the first 
term to scale similarly we require 

KL) = , (29) 

where d is some constant to be determined. Having deduced the required functional 
form of h{L) we now turn to determining the constant d and the precise location and 
properties of the weak peak. 

Taking the first two derivatives of il){n) we have 
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r{n) = - k-— + - . (31) 

Setting il)'{n*) = 0, gives to leading order in L 

n* ^dV(fc+i)^^7(fc+i)(inL)i/('=+^) , (32) 

^in*) - ^^~^^^ lnL-— ^InlnL, (33) 
^ ' k+l k+l ^ ' 

^"{n*) ^ - kL-''/^''+'^ {\nLf-'^^^''+'^ . (34) 

These expressions imply that the height of the peak at n* is 

Pin*) ~ L'^^'i-mk+i) (inL)-''/('+^) , (35) 
and the width of the peak is given by 

An* = \r{n*)\-'^' ~ ^^^^yik^i)/(2ik+i)) ^3^^ 

Thus the weight w of the weak peak in p{n) , which we define as 

w = p{n*)An* (37) 

is given by 

The balance equation ( fT5|) is satisfied by the contribution of the weak peak and 
becomes asymptotically, 

L''-' ~ {n*)''w . (39) 



Inserting the expressions (l38il and ( 1321) implies, ignoring logarithmic factors, 

,.6_l(t±il. (40) 
Thus, from flTSl) . p{n) behaves asymptotically as 

sik + l)\ lnLl'/^'+'^ 



p{n) ~ n exp 



n 



(41) 



k J L } {k + l)L>' 

When s < kh/{k + 1) one has d > as required. Also we require that the weight fl38l) 
of the peak should not diverge as L — 00 otherwise the system would be in the Strong 
High Density Phase. This implies from (138!) and (140|) that s>k/{k + 1). 

To summarise, for k/{k + 1) < s < kh/{k + 1) we have a probability distribution 
which initially decreases from a finite p(0) as a power law as in Critical Phase A, but 
with a weak peak at high n. The weak peak allows the creation-annihilation balance 
condition to be satisfied (see Fig [2]). At this point let us recap the properties of the 
weak peak, ignoring logarithmic correction factors: 

n* ~ An* ~ L^l^^+^) (42) 
pin*) ~ (43) 
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Within the region k/{k + 1) < s < kb/{k + 1) there are in fact two phases 
distinguished by the behaviour of the global density. To see this we note that the 
number of particles per site, Uyjp, associated with the weak peak is 

n^p ~ n*w ~ L^''/i>'+^)-' (45) 

This gives the contribution of the weak peak to the density. Clearly this contribution 
diverges or goes to zero according to whether or not s < 2A;/(A; + 1), thus implying two 
distinct regimes. 

Critical Phase B: 2k/{k + 1) < s < kb/{k + 1) 

In this region n^p approaches zero in the large L limit and the global density of the 
system is controlled by the power law part of p{n). Thus the global density is p = pc 
where pc is the critical density of the corresponding conserving ZRP ([6]). 
Weak High Density Phase: k/ {k + 1) < s < 2k/ {k + 1) 

In this region n^p diverges so the density is dominated by the weak peak and the total 
number of particles increases as ~ ^2fc/(fc+i)-s+i_ 

In these two phases most of the sites of the system form a fluid with typically low 
occupation numbers. In addition a subextensive number of sites, Lw ~ ]^i-s+k/{k+i) 
are highly occupied with n n* particles where n* ~ diverges sublinearly 

with L. We term these highly occupied sites mesocondensates. The expected number 
of such mesocondensates is given by wL ~ L^-s+k/ik+i) ^^^^Yi can vary from 0{L) to 

Q^J^l-k{b-l)/{k + l)y Jf 

6 > 2 + p (46) 

the average number of mesocondensates can be very much less than one. Thus, in this 
case we typically do not expect to observe any mesocondensates. In terms of s, this 
happens when 

kb 2k +1 

> s> . (47) 



k+1 k + 1 

Another distinction between Critical Phase A and the phases characterised by a 
weak peak is the behaviour of arbitrary moments of n 

Pr = J dnn'^p{n) . (48) 

In Critical Phase A pr diverges when r > b — 1. On the other hand, the contribution to 
Pr from a weak peak is {n*Yw ~ which diverges when r > s{k + 1)/A; — 1. 

Therefore in Critical Phase B and the Weak High Density Phase which moments diverge 
depends on the precise values of s,k. This is in contrast to Critical Phase A where 
throughout the phase the same moments {pr for r > b — 1) diverge. 

3.5. Summary of phase diagram 

We have identified the following phases which are illustrated in a typical phase diagram 
in Figured] 
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• k < s — Low Density Phase (LD) 

Here p{n) decays as L~^^~^^''^ and the global density vanishes as L~^^~^\ 

• kb/{k + 1) <s <k — Critical Phase A (CA) 

Here p{n) decays algebraically as with a cut-off at where y = {k — s)/{k — 
6 + 1). The global density throughout this phase is the critical density p^.. 

• 2k/ {k + 1) < s < hk/{k + 1) — Critical Phase B (CB) 

Here p{n) decays algebraically as l/n^ and is cut-off by a weak; peak, with properties 
fH2M). at rf ~ [\Y^Lf/ik+^)L^/i.k+i)_ jj^g global density is the critical density p^. 

• k/{k + 1) < s < 2k/{k + 1)) — Weak High Density Phase (WHD) 

Here p{n) decays algebraically as 1/n'' and is cut-off by a weak peak at n* ~ 
(l^iy/ik+i)];^k/ik+i)^ The global density diverges as p ~ i2k/(k+i)-s 

• s < k/{k + l) — Strong High Density Phase (SHD) 

Here p{n) is sharply peaked around n* ~ and the global density diverges as 

n*. 

Numerical evidence for these phases is presented in Figure [3] and Figure HI In Figure [3] 
typical distributions p{n) for the phases with non-zero density are presented and compare 
favourably with the theoretical predictions. The theoretical predictions were generated 
by taking the predicted asymptotic forms of A and inserting into ( fTOl) . In the case of 
Critical Phase A, the constant g was taken to be 1 although, g could be used as a 
parameter to improve the fit to the simulation data. In Figure HJ p{n) is plotted in 
Critical Phase A for different L illustrating the L dependence of the cut-off. 



3. 6. Nearest Neighbour Hopping 

As noted above we expect the mean-field approximation to be applicable to the fully- 
connected lattice in the limit L — * oo. We also carried out numerical simulations of the 
model on a one- dimensional lattice with totally asymmetric hops to nearest neighbour 
sites and periodic boundary conditions. The results are given in Figure [3] and Figure H] 
and compare well with the mean-field predictions. 

3. 7. The case h<2 

As noted in the introduction, in the conserving model with b < 2 condensation does 
not occur. For the non-conserving model the mean-field analysis follows closely that 
presented above. One finds that Critical Phase B no longer exists and Critical Phase A 
becomes a Weak High Density Phase. Thus for b < 2 one no longer has critical phases 
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Figure 3. Steady-state distributions of the number of particles on a site from 
simulations of the non-conserving ZRP model on a fully connected lattice (open shapes) 
and a Id lattice (character symbols), compared with theoretical curves (dashed lines). 
Simulations were run on a system with L=5000 lattice sites and b = 2.6, fc = 3. Data 
are shown for: Critical Phase A, s = 2 (o, + and - -); Critical Phase B, s = 1.7 (□, x 
and - -); Weak High Density Phase, s = 1.2 (O, * and - -); Strong High Density Phase, 
s ^ 0.4 (A, • and - -). 
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Figure 4. Steady state probability distributions from simulations in Critical Phase 
A, illustrating L dependence of the cut-off. Here b = 2.6, A; = 3 and s ~ 2, with 
L = 1000, 2000, 4000, 8000. The legend indicates for each data set the system size 
and the geometry used: fully-connected lattice (FC) or onc-dimensional lattice (Id). 
The dashed line corresponds to P{n) ^ n^^ . 
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and the phase diagram reduces to Low Density, Weak High Density and Strong High 
Density Phases with phase boundaries given hj s = k and s = k/{k + 1). 

4. Discussion 

In this work we have studied a ZRP with non-conserving dynamics by means of a mean- 
field theory. For the choice of rates dSHSj) we find a rich phase diagram with five distinct 
phases. For low creation rate we find a low density phase with exponentially decaying 
single-site occupation distribution. On the other hand at high particle creation rate we 
find a strong high density phase where the single-site occupation is sharply peaked at a 
large occupation which diverges with system size. 

The most interesting phases are in the intermediate region. Here we find two distinct 
critical phases and a weak high density phase. In the Critical Phase A the occupation 
distribution decays algebraically with occupation number with a finite-size cut-off which 
diverges with system size. On the other hand in Critical Phase B there exists a weak but 
broad peak at large occupation in addition to the algebraic decay. The height and width 
of this peak scale algebraically with system size and the area under the peak vanishes. 
This peak corresponds to a large, but subextensive, number of mesocondensates each 
containing a large but subextensive number of particles. The contribution of this peak 
to the global density vanishes in the thermodynamic limit leaving the global density as 
Pc- In addition to these two phase we find a Weak High Density Phase whose structure 
is very similar to Critical Phase B except that the contribution of the weak peak to the 
global density is dominant and the global density diverges. 

Since the ZRP is a prototypical model for condensation phenomena, we might 
expect the phases established here to be displayed in other driven systems. It would 
be of interest to explore this possibility by studying other microscopic models. For 
example p{n) given in (fTOj) is related to the single-site weight that would be obtained 
in a conserving ZRP with hop rate given by u{n) = u{n) + a{n). This would imply a 
conserving ZRP with non monotonic hop rates. Studies of this conserving model have 
revealed that the phenomenon of multiple condensates exists there as well [20] . 
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